Quantum speed limits in open system dynamics 
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Bounds to the speed of evolution of a quantum system are of fundamental interest in quantum metrology, 
quantum chemical dynamics and quantum computation. We derive a time-energy uncertainty relation for open 
quantum systems undergoing a general, completely positive and trace preserving (CPT) evolution which pro- 
vides a bound to the quantum speed limit. When the evolution is of the Lindblad form, the bound is analogous 
to the Mandelstam-Tamm relation which applies in the unitary case, with the role of the Hamiltonian being 
played by the adjoint of the generator of the dynamical semigroup. The utility of the new bound is exemplified 
in different scenarios, ranging from the estimation of the passage time to the determination of precision limits 
for quantum metrology in the presence of dephasing noise. 
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How fast can a quantum system evolve? Quantum mechan- 
ics acts as a legislative body imposing speed limits to the evo- 
lution of quantum systems. While these limits are both ulti- 
mate and fundamental, at the same time, their existence is at 
the center of a surge of activity, as a result of their manifold 
applications, including the identification of precision bounds 
in quantum metrology [ ], the formulation of computational 
limits of physical systems [ ], and the development of quan- 
tum optimal control algorithms [ ] . 

Bounds on the speed of evolution are intimately related to 
the concept of the passage time i>, which is the required time 
for a given pure state \%) to become orthogonal to itself under 
unitary dynamics [ ]. One of the early answers to this problem 
was provided by Mandelstam and Tamm (MT), who showed 
that the passage time can be lower-bounded by the inverse of 
the variance in the energy of the system so that 
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where AH = ((H 2 ) — {H) 2 ) x l 2 , whenever the dynamics un- 
der study is governed by an Hermitian Hamiltonian H [5-8]. 
A simple geometric interpretation of this result was provided 
by Brody using the Fubini-Study metric in the Hilbert space 
spanned by the initial state and its orthogonal complement [ ] . 
Indeed, the passage time problem can be posed as a quantum 
brachistochrone problem. From this perspective, a particu- 
larly exciting result was found: whenever the Hamiltonian is 
non-Hermitian PT-symmetric, the passage time can be made 
arbitrarily small without violating the time-energy uncertainty 
principle [10, 11]. A second bound, due to Margolus and Lev- 
itin (ML), takes the simpler form T > | (h)-^ wnere me zero 
of energy is generally shifted to the ground state energy so 
that Eq =0 [12]. This bound has been applied to ascertain 
fundamental computational limits in nature [2, 13]. 

Despite the growing body of literature on the subject, the 
analysis has almost exclusively been focused on unitary dy- 
namics of isolated quantum systems. An analogous bound for 



open quantum systems is highly desirable, since ultimately 
all systems are coupled to an environment [14, 15]. As an 
example, such a bound on the evolution of an open system 
would help to address the robustness of quantum simulators 
and computers against decoherence [16], as well as the rel- 
evance of the specific nature of the noise, and in particular 
whether or not it is Markovian, in phase estimation problems 
of interest in metrology and precision spectroscopy [17, 18]. 

The MT bound can be derived by considering the time evo- 
lution of the overlap a = \(Wt\¥o)\ between the initial state 
|y/b) and the quantum state \ y/ r ) at time t subject to a unitary 
evolution U (t) — exp{— iHt/h}. It can be shown that the MT- 
limit (eq. 1) is achievable, as for a suitable Hamiltonian H we 
can satisfy the differential equation fi^r = —ZAHaV 1 — a 2 
which for a = cos (j) is easily seen to result in h(j) = AH thus 
matching the MT bound [19]. 

In the case of open system dynamics we need to consider 
general non-unitary quantum evolutions and have the free- 
dom to choose a variety of distance measures between quan- 
tum states. One natural choice here is the fidelity between 
two mixed states p and a, which is given by F(p,a) = 
yfp Oyfp]. The quantum speed limit then provides a 
lower bound on the time T that is required to achieve, for a 
given initial state p(0) and a target fidelity f targe t, the condi- 
tion F(p,,po) < /target subject to an open system evolution. 
Ideally, such bounds should reduce to the MT bound in the 
case of unitary dynamics on pure states and/or be easy to com- 
pute. 

Bounds on T may be derived by taking inspiration from 
the variational characterization of the fidelity F(p s ,a ) = 
max[|(v/ S£ |0 S£ )|] [20], where the maximization is over all 
I y SE ) (\<I> SE )) on a larger Hilbert space 3? SE that are purifica- 
tions of the mixed states p s (a s ) on the smaller system S, that 
istr £ [|y/ S£ )(y/ S£ |] =P S (tr £ [|0 S£ )(0 S£ |] = ff s ). Then for any 
specific purification the inequality F(p s ,o s ) > \ (\ir SE \<j) SE )\ 
holds. A general time evolution of a subsystem p s = ££p s 
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can always be generated by a joint unitary dynamics U t SE of 
the system with an environment such that 
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where p^ = tr B [| Yq E ) (Vo^ll- Such a dynamics will be gener- 
ated by a suitable Hamiltonian Hf E but it should be noted that 
the choice of | y/g £ ), U t SE and thus Hf E is not unique. 
Now we can make use of the fact that 

ftarg* > HPt'Po) > F(U t SE \xi4 E ),\xi4 E )) 

for any choice of purification of pf and Pq and any choice 
of unitary dynamics U t SE that generates p s = J>fp s on the 
subsystem. This implies that any choice of purification and 
unitary evolution will achieve f target > F(U t SE \ V/q £ ), | Wq E )) 
at an earlier time t than f tur get > ^(Pt>Pq)> f < T - As a 
consequence, for any choice of |i//q }, £/, and thus // r S£ we 
obtain a lower bound on T. If Hf E is given, then we can com- 
pute AHf E in the state |y/f £ ) and immediately provide a lower 
bound on x via the MT and M0 = hf^ds^ < J^dsAH^. 
Needless to say, performing the optimization over all possible 
purifications and all possible Hf E is a challenging task that 
will be very hard to perform in the general case. Two routes 
are suggested themselves. Firstly, well chosen |W? ), Uf E and 
thus Hf E will lead to excellent bounds for reasonably simple 
cases. Secondly, analytical lower bounds on x may also be ob- 
tained by studying different distance measures that are easier 
to handle and thus admit closed formulae for lower bounds. 

Here we follow this second approach to find an analytical 
and easy to compute lower bound on the speed of evolution in 
open quantum systems. We shall derive a bound analogous to 
the seminal result by MT where the energy variance of the ini- 
tial state is replaced by a more general measure taking into ac- 
count the coupling to the environment. We shall pay particular 
attention to the dynamics governed by a dynamical semigroup 
in which case the evolution of the system is ruled by a mas- 
ter equation of the Lindblad form [ ] . It has recently been 
pointed out that the MT and ML bounds are violated under 
non-unitary dynamics [22]. We shall nonetheless show in the 
following that Markovian systems are subjected to a MT-type 
of bound where the adjoint of the generator of the dynami- 
cal semigroup plays the role of the system Hamiltonian in the 
unitary case. 

Decay of an open quantum system Consider a given sys- 
tem described by a state po (from now we drop the upper in- 
dex S for convenience) coupled to an environment in a state 
p E , and assume both system and environment are weakly cou- 
pled such that the initial global state can be approximated 
by po (8 Pq- Let the global reversible dynamics be governed 
by a unitary evolution operator U,. The reduced dynamics 
of the system is given by a one-parameter family of dynam- 
ical maps p H> Y r p := \ie[U,Pq& p E u}], parameterized by 
the time variable t E M. + . Whenever the typical time scale 
of the environment is much smaller than that of the system, 
one can assume a Markovian dynamics. Under Markovian 
dynamics, such maps form a quantum dynamical semigroup 



ft+sP — "ft%P> t,s > (we assume that the open system is 
not subjected to an external time-dependent field so that the 
generator of the quantum dynamical semigroup is time inde- 
pendent). Any such map can be represented by a Markovian 
master equation 



dp, 
dt 



where the generator of V, admits the Lindblad form [ ] 



J?p = --[H,p] 
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such that YtPo — e'^Po- In such scenario we might pose the 
following question: Which is the bound to the speed of evo- 
lution from an initial state po under the action of a quantum 
dynamical semigroup X? To answer this question we intro- 
duce as a figure of merit the so called relative purity [ ] 



fit) = ^ 
t{) tr(p2) 



(4) 



which is a generalization of the survival probability y{t) = 
\(x\e~ iH,/t '\x)\ 2 often used for a pure state \%) subject to a 
Hamiltonian H, and that has proved useful in studying quan- 
tum speed limits in the unitary case [13]. 

Derivation of the bound from the (Lindblad) master equa- 
tion Let us now characterize the decay rate of the relative 
purity. Note that whenever the generator admits a Lindblad 
form (i.e. for a Markovian quantum master equation), 
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tt|poJS?p,] _tr[-S?W f ] 



(5) 



tr(p 2 ) tr(p2) 
where the adjoint of the generator of the dynamical map reads 

if f Po = j. [H, po] + £FtpoF* - \ {FkFhPo } ■ (6) 

k 

The rate of change of / can then be bounded using 
the Cauchy-Schwarz inequality for operators, |tr(A'B) 
tr(A t A)tr( J S t J S). Then 



< 



< VM^Po^Mp/l/trpo 2 
v /tr[(^tp )2 ]/trp 2 ; 
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that is, by making reference exclusively to the initial state and 
the dynamical map. Let us parametrize f(t) — cos# with 
& € [0,7r/2], Upon integration between # = (/(0) = 1) 
and a final & = 9, the following bound to the required time of 
evolution is found 



Icos0 - lltrp,? 
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This generalizes the MT uncertainty relation for open quan- 
tum systems governed by a Markovian quantum master equa- 
tion. The generalization to a time-dependent Lindbladian 
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_£f (t) is straightforward and reads 



T 6 > 



(9) 



-- 1 r^xdt. 



where X — Tg 1 J 

Derivation of the bound using general quantum channels 
To remove the Markovian approximation, we note that any 
kind of time evolution of a quantum state po can be written 
in the form p, = Y,a K a{t^)poKa{t,Q). In particular, K a (t,Q) 
is independent of po if the dynamical map is induced from an 
extended system with the initial condition p^ E = po <E> p £ (0). 
Then, the dynamical map is said to be universal. Let such map 
govern the evolution and consider 

f(t) =ti\p oPt ] = J>[p OJ M?,0)po4(r,0)]. do) 

a 

Parametrizing /(f) = cos 6, a bound can be derived 
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where ||A|| = ^/tr(A"A) is the the the Hilbeit-Schmidt norm 
of A. Details of the derivation are provided as an Appendix. 

Applications The bound to the speed of evolution pre- 
sented above is the main result of this paper. In the following 
we shall analyze some particular cases to illustrate its use: 

Passage time. - Under unitary time evolution, the passage 
time is the minimum time required for a time evolving state 
\x(t)) to become orthogonal to its initial value |j(0)). Let 
us consider a pure state such that trpg = 1 and let J>f^ p = 
-J^p = i[H, \x) (X\] A ^ follows from Eq. (9), that 



1Q > 



h 
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Alternatively, for a = l,K = exp[-i(H- (x\H\x))t /h], x e > 
a factor 1 /\/2 smaller. The usual definition of the pas- 
sage time %p = 2^=, refers to the orthogonalization measured 
by the square root of the fidelity (sometimes referred to as the 
integrity or survival probability amplitude), which for a pure 
states under unitary dynamics reduces to |(^(0)|^(f))| [4-7]. 

Non-Hermitian Hamiltonians. - Non-Hermitian Hamiltoni- 
ans are ubiquitous in quantum physics and enjoy of a wide 
range of applications from quantum optics [24] to reactive 
scattering [ ]. Their standard derivation is based on Fesh- 
bach's partitioning theory, that allows to describe the effective 
dynamics of a quantum system governed by a Hamiltonian H, 
when restricted to a given subspace associated with projector 
P (with complement Q, such that P + Q = 1 , P 2 = P, Q 2 = Q). 
The effective Hamiltonian governing the dynamics in the re- 
stricted subspace, H efl = PHP + PHQ(E - QHQ) l QHP, 
is generally non-Hermitian. Under Hp the density matrix 
ip = (H e ffP — pH^ s )/h. Similarly, in open systems under 
Markovian dynamics it is customary to split the generator 



of the dynamical map in two contributions «5f £ and 3>, i.e. 
«5f = JS^ + £F. J£ c describes the coherent evolution associated 
with the non-Hermitian Hamiltonian H e ft = H — ih\ L^/ft, 
while the dissipator Sp = Y^k^kPF^ is associated with spon- 
taneous decay, and it is a jump operator [ ]. More generally, 
let Heff = H — ;T, where H and F are both Hermitian opera- 
tors, so that Jzf t p = —i[H,p]/h — {r,p}//j. Noting that upon 
setting S>p = 0, the bound to the speed of evolution under 
non-Hermitian Hamiltonians still holds, it follows from Eq. 
(7) that 
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where the last line applies exclusively to pure states. Using 
Eq. (11) with a = l, K = exp[-i(H-(x\H\x))t/h-(T- 
i(x\T\x))t/h], one finds 1 j\fl times the same expression. 

From quantum speed limits to metrological bounds The 
ultimate bound to parameter estimation is dictated by the abil- 
ity to efficiently discriminate neighboring quantum states. In 
a seminal paper [ ], Braunstein and Caves (BC) derived a 
quantum Cramer-Rao bound for the uncertainty in the (local) 
estimation of a classical parameter (j) of the form: 



A<j> > 



1 



(14) 



where Fq denotes the quantum Fisher information and V is 
the total number of repetitions of the experiment where a 
(j) -dependence is linearly imprinted via a general evolution. 
When the dynamics is unitary, an initial preparation of a probe 
state in a cat (GHZ) state of Af subsystems allows to satu- 
rate the lower bound and achieve a Heisenberg-limited reso- 
lution where A<p ~ 1 /N. If the subsystems are used inde- 
pendently, so that the input state is factorizable as product 
states, only the standard scaling dictated by the central limit 
theorem A0 ~ 1 /y/N is achievable. This implies that the er- 
ror bars in the actual estimation of a parameter (j) could be 
reduced by 1 v^V by means of employing an entangled input 
probe provided that the system evolves unitarily. Whether or 
not the standard scaling can be surpassed when the system's 
dynamics dynamics is open is a most relevant issue where 
only partial results are known. Motivated by experiments on 
precision spectroscopy, where a phase difference is estimated 
which is proportional to the detuning between an external os- 
cillator and a selected atomic frequency, we will focus here on 
phase estimation problems under dephasing noise. Assuming 
decoherence to be Markovian and affecting each subsystem 
independently (local noise assumption), it was shown in [ ] 
that this type of noise renders product and maximally entan- 
gled states metrologically equivalent, and argued that Marko- 
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vian dephasing would restore the standard scaling with an op- 
timal resolution to be achieved by a type of partially entangled 
states so that At))" 1 * /Aty 1 ' = \/\fe. Subsequent work proved 
this bound to be achievable asymptotically [ ] but only very 
recently it was proved in all generality that the bound is sharp 
and coincides with the one imposed by the maximization of 
the quantum Fisher information [ ]. The metrological equiv- 
alence of product and maximally entangled state preparations 
under Markovian decoherence can be predicted with the new 
bound eq. (8), which yields the ratio ?ghz = t p /N, where Jqhz 
and t p are the optimal interrogation times when using maxi- 
mally entangled and product state inputs, respectively. This 
can be easily shown by writing the dephasing master equation 
in the interaction picture as 

p = -YP + Yo z po z , (15) 

and considering a pure state po = |#)(#| of the form \%) = 
(|0> + |1»/V2. Then, 

^Po = -ypo + yc; z p<)C> z (16) 

and 

ti[(^p ) 2 ]=2 Y 2 . (17) 

This yields a minimal orthogonalization time t p = l/y/ly. 
Repeating the same procedure for a maximally entangled in- 
put of the (GHZ) form p = \x)(x\ with \%) = (|0}® iv + 
\\)® N )/\/2, we obtain an optimal interrogation time ?ghz = 
t p /N which leads to A0 GHZ = A(j) p when the resolution is es- 
timated operationally as A(j> =< AO 2 > 2 /*fv | d<°> |, with 
O denoting a projective population measurement, which is 
known to be optimal for this specific context. Alternatively, 
we can estimate the Fisher information in the form 

where pi = tr(p^,P ( ) and P, is a population projective measure- 
ment. Note that this measurement procedure is optimal in this 
context. The resulting expressions for product and cat states 
are, respectively 

F p = Ne- 2rt t 2 , (19) 
Fghz - N 2 e- 2N y t 2 . (20) 

The ratio A0 GHZ /AtyP = ^(v p F p )/(v ghz F Gii z) therefore 
equals 1 when considering the optimal interrogations times as 
dictated by the bound eq. (8). Moreover, for pure states po = 
|^)(^| and the case of Markovian pure dephasing Jzf^Po = 
rE thPo + ofWi* j] we have that ^[(j^p p) 2 ] = 

yyj (Ni-2NZ k \(xWP\X)\ 2 +L kl M°P°P\x)\ 2 < 
^/2yN (Note that this may be generalized to the mixed state 
case and any form of local noise as the locality implies that 
number of terms in Jzf^po grows linearly in the number of 



subsystems N). Then with eq. (8) and the fact that the Fisher 
information obeys F < N 2 [29, 30], the limit on the speed 
of evolution imposes the persistence of the standard scaling 
A(j> ~ 1 /y/N no matter how weak the dephasing rate. This is 
a result that is now firmly established [28, 31] and that comes 
out in a rather natural fashion within this new framework. 

So far we have exploited specifically the fact that the sys- 
tem's dynamics is ruled by a Lindblad master equation. How- 
ever, our initial derivation considered a dynamical map that is 
CPT but not necessarily divisible. As a result, the bound could 
be valid for non Markovian dynamics that admit a representa- 
tion in terms of a CP map. We have evaluated the prediction 
for the optimal interrogation times of product and cat states 
for a model of non-Markovian dephasing of this type, as pro- 
posed in [32], and obtained the ratio fGHZ = tp/N, just as in 
the Markov case. This seems to be in contradiction with re- 
cent results for models of non Markovian dephasing, which 
predict a ratio Jghz = t p / VN and raises an interesting conjec- 
ture with which we finish this section. There could exist forms 
of coloured noise for which the metrological equivalence be- 
tween cats-products input probes still holds. This inequiva- 
lence in the achievable resolution of a phase estimation could 
then be exploited to quantitatively quantify non-Markovianity. 

Conclusions- A bound to the speed of evolution under 
an open-system dynamics has been provided, generalising the 
classic result by Mandelstam and Tamm known for the unitary 
case. In the Markovian limit, we have shown that the adjoint 
of the generator of the dynamical semigroup plays the role of 
the commutator with the Hamiltonian in the MT bound. De- 
spite the fact that the bound is not tight, in the sense of non 
coinciding with the unitary solution for closed systems, it al- 
lows to naturally predict the inaccessibility of the Heisenberg 
limit under Markovian noise. Moreover, when using the gen- 
eral form of the bound for universal channels, the new limit 
on the speed of evolution suggest the inequivalence of differ- 
ent forms of coloured noise for precision spectroscopy. Our 
results are applicable to a wide variety of scenarios including 
bounding decoherence rates [33], and quantum speed limits in 
dissipative state preparation [34], quantum computation and 
simulation assisted by dissipation [35]. 
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Derivation of the bound using quantum channels 

Any kind of time evolution of a quantum state po can be 
written in the form p t = L a ^«(f ,0)po^a(f ,0). In particular, 



K a (t,Q) is independent of po if the dynamical map is induced 
from an extended system with the initial contidion p^ £ = po <S> 
Pq. Then, dynamical map is said to be universal. Let such 
map govern the evolution and consider 



m = ^ 
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£tr[p OJ Mf,0)po4(r,0)]. (2D 



trpo tr Po a 
For compactness, let us denote K a = K a (t,0). It follows that 



t{ ) " trpg ' 



£ tr[p (4po4 + K aPo Kl] . (22) 



Taking the absolute value at both sides, 



1/(0 1 < -jL= £ [|tr(p 4po4)l + |tr(p ^po4) I] • 
Mpo 2 ] a 
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Using the Cauchy-Schwarz inequality |tr(AB)| < 
[tr(A t A)tr(B t B)]i, 



1/(01 < 



£(A/tr[4po4^aPo4 



tr[Kap Q k r a k a poKl,}), 
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We next use the fact that | / dtf(t) \ < J dt\f(t)\. Parametriz- 
ing f(t) = cos 9, a bound can be derived 



T > 
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where X = T Q 1 J^Xdt and ||A|| = ^/tr(AtA) is the Hilbert- 
Schmidt norm of A. In the second line, we have used | cos 9 — 
1| >49 2 /7t 2 . 



